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BOUNDARY  LAYER  FLOW  WITH  HEAT  TRANSFER 
AT  SUPERSONIC  AND  HYPERSONIC  SPEEDS 


Prepared  by: 
Jerome  Persh 


ABSTRACT:  A  theoretical  investigation  of  compressible  turbulent 
boundary  layer  flow  with  and  without  steady  state  heat  transfer 
has  been  conducted.  This  investigation  is  based  on  a  ->  •  .;plc 
physical  model  of  the  flow  suggested  first  by  Prandtl  and  used 
later  by  Donaldson.  The  physical  model  consists  o<  a  laminar 
sublayer  region  with  a  linear  velocity  profile  and  an  outer  tur¬ 
bulent  portion  with  a  power  law  velocity  profile.  Comparisons 
between  theory  and  experiment  demonstrate  that  the  analysis 
yields  good  results  for  compressible  turbulent  boundary  layer 
flow  with  and  without  steady  state  heat  transfer. 
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This  report  contains  the  results  of  a  theoretical  investi¬ 
gation  of  compressible  turbulent  boundary  layer  flow  with  and 
without  steady  state  heat  transfer.  The  significance  of  this 
work  is  apparent  when  it  is  considered  that  although  a  great 
deal  of  experimental  and  theoretical  information  exists  for 
supersonic  turbulent  boundary  layers  in  the  absence  of  heat 
transfer,  there  are  relatively  few  detailed  investigations 
in  the  supersonic  and  hypersonic  speed  ranges  that  include 
the  effects  of  heat  transfer. 
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SYMBOLS 

Cf  «  Local  skin  friction  coefficient  based  on  tree-stream 
conditions,  2  t'w/^*00 2 

Cf  «  Mean  skin  friction  coefficient  based  on  tree-stream 
conditions,  2D 

f'n'xi  * 

cji  <•'  Incompressible  local  skin  friction  coefficient  for 
zero  heat  transfer  based  on  free-stream  conditions 

Cpi  *  lncospressible  sean  skin  friction  coefficient  for  zero 
heat  transfer  based  on  tree-stream  conditions 

D  -  Drag  force 

H  *  Boundary  layer  shape  parameter,  sV© 

k  -  Constant  in  nixing  length  law 

j£  -  Mixing  length 

M  -  Mach  nunber 

n  -  Exponent  in  power  law  velocity  profile  representation 

r  •  Ratio  of  total  shear  stress  to  viscr’;a  shear  stress 

r.f.  «  Recovery  factor 

Re  -  Reynolds  nunber 

T  -  Local  static  temperature 

u  «■  Mean  velocity  component  in  x-direction 

u*  -  Velocity  parameter,  u/ur  (based  on  wall  conditions) 

ur  -  Friction  velocity 

x  “  Axial  distance  along  surface 

y  -  Distance  perpendicular  to  surface 

y+  -  Wall  distance  parameter,  iur/j  (based  on  wall  condi¬ 
tions)  1 
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x  -  Values  based  on  distance  from  leading  edge  of  plate 

S  »  Va.1  s  based  on  boundary  layer  thickness 

8  -  Values  based  on  boundary  layer  momentum  thickness 

-  Values  based  on  free-stream  conditions  outside 
boundary  layer 
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A  THEORETICAL  INVESTIGATION  OF  TURBULENT 
BOUNDARY  LAYER  FLOW  WITH  HEAT  TRANSFER 
AT  SUPERSONIC  AND  HYPERSONIC  SPEEDS 


INTRODUCTION 

1.  Despite  a  lack  of  experimental  data,  numerous  formulae 
have  teen  developed  for  the  variation  of  turbulent  skin 
friction  on  a  flat  plate,  with  and  without  steady  state 
heat  transfer.  The  reports  of  Rubesin,  Maydow,  and  Varga 
(reference  a),  and  Chapman  and  Kester  (reference  1;)  include 
good  resumes  of  several  theoretical  treatments  of  this  proi.lem. 
All  of  the  analyses  reviewed  in  these  references  make  use 

of  empirical  constants  which  are  drawn  from  incompressible 
experimental  data.  Recent  experimental  results  (reference  c) 
have  demonstrated  that  the  empirical  incompressible  constants 
utilized  are  affected  by  heat  transfer.  Specifically,  it 
has  been  found  that  the  assumption  that  the  edge  of  the 
laminar  sublayer  occurs  at  fixed  values  of  the  parameters 
u£  (or  y£)  is  not  strictly  valid.  Experimental  results 
indicate  that  the  value  of  u£  does  not  only  vary  with  heat 
transfer,  but  to  some  extent  with  Reynolds  number  and  Mach 
number.  It  was  felt,  therefore,  that  a  theoretical 
approach  which  is  based  on  a  realistic  physical  model  of  the 
flow,  and  which  allows  a  prediction  of  the  quantities  at  the 
edge  of  the  laminar  sublayer,  is  expedient  at  the  present 
time , 

2.  Such  an  approach  was  originally  devised  by  Prandtl 
(reference  d)  and  recently  extended  to  compressible  flows 

by  Donaldson  (reference  e)  The  physical  model  of  a  turbulent 
boundary  layer  proposed  by  .iiose  Investigators  may  be  briefly 
described  as  follows:  It  is  assumed  that  the  turbulent 
boundary  layer  velocity  profile  can  be  divided  into  two 
regions;  the  wall  adjacent  region  called  the  laminar  sublayer, 
whore  the  velocity  varies  linearly  with  distance  from  the 
surface,  and  the  outer  turbulent  portion,  which  is  represented 
by  a  power  profile.  The  intersection  of  these  two  profiles  is 
defined  as  the  edge  of  the  laminar  sublayer. 

3.  It  is  the  purpose  of  this  investigation  to  oxtC'tid  and 
revise  the  analysis  of  Donaldson  in  order  lo  o.  Lain  consistency 
with  the  most  recent  ami  reliable  experimental  results  for 

low  speed  turbulent  boundary  layers.  The  applicability  oi  this 
analysis  for  compressible  turbulent  boundary  layers  with  and 
without  steady  state  heat  transfer  is  demonstrated  . y  comparisons 
with  supersonic  and  hypersonic  experimental  results. 
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ANALYSIS 

J r..ompressible  Turbulent  Boundary  Layers 

4.  It  hr.<»  been  established  by  numerous  investigators  that  the 
velocity  profile  in  the  outer  turbulent  portion  of  the  boundary 
layer  may  bo  adequately  represented  by  a  power  profile  of  the 
form  i 


while  experimental  evidence  indicates  that,  in  the  laminar  sub¬ 
layer  region,  the  velocity  profile  is  essentially  a  straight 
line 

H_  *>  1 —  x  constant  o  —  y  —  S 

UL  6  L  L 


The  boundary  layer  is  thus  divided  into  a  turbulent  portion 
described  by  Eq.  (1)  and  a  laminar  region  having  a  linear  ve~ 
locity  profile  (Eq.  2).  This  is  shown  in  Figure  1  with  the 
real  conditions  in  the  transition  region  indicated  by  a  dashed 
line. 

5.  Several  general  relations  for  the  local  skin  friction  coef¬ 
ficient  can  be  deduced  using  the  preceding  postulates  regarding 
the  boundary  layer  velocity  profile  together  with  various  as¬ 
sumptions  regarding  the  shear  stress  at  the  edge  of  the  laminar 
sublayer.  These  relations  necessarily  embody  unknown  functions 
which  must  be  evaluated  empirically. 


6,  Donaldson  (reference  e)  introduced  an  empirical  constant 
relating  the  total  shear  stress  and  the  laminar  shear  stress 
in  order  to  compute  the  skin  friction.  Taking 


T. 

Lam . 


I 


r, 


Turk . 


t  -  constant 


(3) 


and  evaluating  r  at  y  *  from  the  power  profile  given  by 
Eq.  (1),  Donaldson  (reference  e)  derived  the  following  rela¬ 
tion  for  the  skin  friction  coefficient 


2 


1) 


***** 
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and  evaluated  the  constant  (r-l)A2  empirically  using  Rlasius  ’ 
(reference  d)  skin  friction  lav.  In  Donaldson's  analysis,  the 
velocity  profile  exponent  (n)  is  considered  constant,  and  the 
Prandtl  nixing  length  lav 

(t  -  k  y  (5) 

is  arsumed  in  order  to  calculate  the  turbulent  shear  stress, 
rTurb.* 

7.  In  the  present  report  the  velocity  profile  exponent  (n)  vas 
taken  as  a  variable  and  the  method  by  vhich  Eq.  (4)  vas  correlated 
vith  experiment  is  es  follows:  First,  Eq.  (4)  vas  equated  to  the 
Karman-ftchoenherr  incompressible  skin  friction  law 

c.jjL  **  0,0568 _ 

[Logic  <2Ree>]  t^io  (2Ree)  +  °-868]  <6) 


vhich  is  regarded  as  a  good  representation  of  incompressible 
turbulent  boundary  layer  skin  friction  coefficients  over  a  vide 
Reynolds  number  range.  This  procedure  yielded  a  variation  of 
n  with  Reg.  A  comparison  vas  then  made  between  experimental  data 
and  the  deduced  variation  of  n  with  Reg.  This  approach  seemed 
logical  becau-e  the  experimental  variation  of  n  with  Reg  is  veil 


\ 


♦  An  examination  of  Figure  2  suggests  that  since  the  shear 
stress  is  nearly  constant  near  the  wall,  an  alternate  assump¬ 
tion  regarding  the  relationship  between  the  laminar  and  turbulent 
shear  stress  may  be  made.  Within  the  sublayer  the  laminar  she.? 
stress  predominates  while  outside  the  sublayer  the  turbulent 
shear  stress  predominates,  and  since  the  transition  region  is 
neglected  in  the  model  for  the  velocity  profile,  the  values  of 
^*Turb.  calculated  from  Eq.  (1)  and  calculate^  fror;  Eq.  (2) 

may  be  taken  as  each  equal  to  the  total "sheer  stress  at  y  -  6,  . 
Thus  a  logical  relationship  between  T  jurb,  and  ^Lam.  would  ' 
appear  to  be 

rTurb.  "  '^La».  y  "  6L 

if  ^Turb.  ^Lam.  *re  computed  as  indicated  above.  In 

addition,  if  the  von  Karman  mixing  length  formula 


£ 


k 


du 

2£_ 

d2u 


is  used  instead  of  the  Prandtl  mixing  length  formula,  each  o» 
the  shear  stress  assumptions  will  lead  to  another'  skin  f’lction 
relation. 


A  study  was  therefore  made  to  determine  which  of  the  as¬ 
sumptions  and  mixing  length  formula  yielded  a  skin  friction 
law  which  gave  the  best  overall  agreement  between  theory  and 
experiment.  It  was  found  that  Eq.  (4)  given  by  Donaldson, 
could  best  be  adapted  to  the  experimental  results.  The  dotnils 
of  this  analysis  are  contained  in  Appendix  A. 
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described  by  a  large  amount  of  experimental  velocity  profile 
data.  It  was  found  that  using  this  procedure  and  taking  £  "  * 

k4 

as  constant  resulted  in  a  variation  of  n  with  Re-  which  is 
a  good  average  of  the  available  experimental  data  (references 
c,  f,  g,  and  h).  This  is  shown  in  Figure  3.  The  value  of  the 
constant  r  -  1  compatible  with  the  n  variation  with  Re.  shown 
~T c2-  * 

in  Figure  3  is  20.0.  This  is  only  slightly  different  from  the 
value  of  22.5  given  by  Donaldson  (reference  e).  Since  in  his 
analysis  the  n  variation  with  Reynolds  number  was  not  considered 
and  the  Blaslus  (reference  d)  skin  friction  law  was  used  to  ob¬ 
tain  the  constant,  it  is  to  be  expected  that  a  slightly  different 
constant  would  be  obtained.  The  incompressible  portion  of  the 
present  analysis  therefore  represents  an  extension  of  the  Donaldson 
analysis  in  that  the  n  variation  with  Reynolds  number  is  con¬ 
sidered. 


Compressible  Turbulent  Boundary  Layers 

8.  Since  for  compressible  flow,  the  temperature  varies  across 
the  boundary  layer,  it  is  apparent  that  the  assumption  of  con¬ 
stant  shear  stress  through  the  laminar  sublayer  violates  the 
stipulation  that  the  velocity  varies  linearly  with  distance 
from  the  surface.  It  is  assumed  that  this  incompatibility  does 
not  introduce  any  serious  errors  in  the  skin  frictica-  results. 
The  subsequent  comparisons  between  theory  and  experiment  tend 
to  confirm  this  assumption. 


9.  The  extension  of  the  foregoing  analysis  to  compressible 
flows  is  straight Torward  and  the  equations  take  the  same  fora 
as  those  given  in  reference (e) .  If  the  value  of  the  constant  r 
obtained  by  the  procedure  described  above  is  assumed  to  be  the 
same  for  both  incompressible  and  compressible  flows  with  and 
without  heat  transfer,  then  Eq.  (4)  is  valid  for  these  cases  if 
the  density  and  viscosity  contained  in  the  Reynolds  number  are 
evaluated  at  the  edge  of  the  lanjinar  sublayer.  It  is  desirable, 
however,  that  the  Reynolds  number  be  expressed  in  terms  of  free- 
stream  properties.  This  transformation  is  presented  in  reference 
(e)and  in  Appendix  B  where  it  is  shown  that  the  resulting  equa¬ 
tions  necessary  to  determine  local  skin  friction  coefficients 
are  as  follows: 


cf 


2 


20  n 


1  -  n 

F'T~n 


L 


-  2  (J-  1 
n  +  1 


(7) 


Too 


1+  r.f  , 


t-  1 
2 


tfj 


Tw  ~  T« 


oo 


(8) 
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10.  Values  of  n  for  use  in  these  equations  are  obtained  from 
the  curve  of  Figure  3.  Using  this  single  curve  as  a  determining 
factor  for  n  implies  that  n  is  uniquely  related  to  Re»  regard¬ 
less  of  Mach  number  or  heat  transfer  condition.  Although  the 
results  shown  in  Figure  3  tend  to  justify  this  postulate,  probably 
the  Mach  number  and  heat  transfer  influences  on  n  are  concealed 
by  the  insensitivity  inherent  in  this  coordinate  system.  A  com¬ 
parison  between  a  plot  of  n  versus  Reg  and  n  versus  Res indicates 
that  some  influence  of  Mach  number  ana  heat  transfer  j.s  probably 
incorporated  in  both  figures,  but  is  appreciably  less  in  the  Reg 
plot.  It  is  felt,  therefore,  that  a  more  accurate  determination 
of  n  may  be  obtained  when  Reg  is  used  as  the  correlating  factor. 


11.  Since  Reg  was  chosen  as  the  correlating  factor  for  n  and 
Eq.  (7)  requires  the  use  of  Re8 ,  a  means  is  therefore  needed  for 
converting  given  values  of  Reg  to  the  Equivalent  Res.  The  needed 
6/6  values  have  been  calculated  using  the  definition  for  6  witn 
Eq.  (1)  and  the  Orocco  temperature  distribution  (reference  e)  for 
a  series  of  Mach  numbers  up  to  20,  a  wide  range  of  heat  transfer 
conditions,  and  n  values  of  5,  7,  9,  and  11.  These  are  tabulated 
in  Table  I.  Also  tabulated  for  the  same  range  of  variables  are 
the  values  of  6*/6  and  H.  It  should  be  noted  that  the  foregoing 
procedure  for  calculating  these  parameters  ignores  the  laminar 
sublayer  because  the  power  profile  is  assumed  to  exist  to  the 
wall.  While  this  procedure  is  not  quite  exact,  it  is  felt  that 
only  small  errors  will  result  because  by  far  the  largest  contri¬ 
butions  to  the  integrals  for  ft*  and  6  occur  outside  the  laminar 
sublayer . 

12.  Using  the  present  analysis,  skin  friction  coefficients  can 
be  calculated  if  the  Reynolds  number  is  given  in  terms  of  either 
the  total  boundary  layer  thickness  (Re  5  )  or  the  boundary  layer 
momentum  thickness  (Rep),  or,  as  will  be  shown  later,  in  terms 
of  the  distance  from  tne  leading  edge.  Since  the  dependence  of 
6/6  with  Mach  number  and  heat  transfer  is  not  considcied  in 
Donaldson's  (reference  e)  analysis,  the  skin  friction  coefficient 
can  be  evaluated  only  if  the  value  of  Re4  is  known. 

13.  Whether  or  not  the  postulated  uniqueness  of  n  with  ile& 
leads  to  serious  errors  may  be  checked  by  comparing  the  influence 
of  n  on  Cf  over  a  range  of  Mach  numbers  for  a  fixed  value  of 

Reg.  A  value  of  Reg  of  8000  was  selected  for  the  check  procedure, 
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and  the  n  value  was  vr.rlcd  between  5.5  and  7.5  which  encompasses 
the  scatter  of  the  experimental  data  at  this  point.  It  was 
found  that  tho  values  of  Cf  over  a  range  of  Mach  numbers  up  to 
10  were  no  more  than  6  percent  above  cv  below  the  curve  drawn 
for  the  value  of  a  theoretically  associated  with  an  Reg  of  8000. 

On  the  basis  of  this  check  it  was  concluded  that  the  skin  fric¬ 
tion  values  obtained  from  the  present  analysis  are  not  particularly 
sensitive  to  the  value  of  n  aacaciAteU  with  the  Reg  In  question. 


14.  This  result  enables  the  approximate  calculation  of  the  varia¬ 
tion  of  cf/cfi  ss  a  function  of  Rex.  If  e.f/cfi  is  not  particularly 
sensitive  to  n,  it  aay  be  assumed  that  for  a  given  Mach  number 
and  heat  transfer  rate,  the  value  of  6/6  along  a  plate  Is  constant. 
This  assumption  is  necessary  to  perform  the  integration  indicated 
in  Appendix  C.  Using  this  assumption,  a  relation  between  Re,  and 
Rex  aay  be  deduced.  The  details  of  this  dsrivation  are  given  in 
Appendix  C,  where  it  is  shown  that  the  resulting  equation  in 


R«»j  -  Rjf*x 


n  +  1 
h  +  3 


X 


U  +  1 

n  +  1 

A  /"  + 

3  \ 

(JL) 

J  1  n  -b 

t rl 

\  20n  / 

n  - 


-  2.52 


Itl 


j 


(10) 


It  is  not  Stipulated,  however,  that  the  value  of  (n)  for  use 
in  equation  (10)  is  a  constant.  Curves  showing  thw  variation 
of  cf/cfi  as  a  function  of  Mach  number  for  several  cor.<*tant 
values  of  wall  temperature  ratio  and  a  constant  value  of  Rex 
of  10?  (Figure  15)  calculated  using  the  equations  of  Appendix 
C  are  in  good  agreement  with  the  empirical  curves  of  Seif' 
(reference  a) .  Results  obtained  using  the  equations  of  Appen¬ 
dix  C  should  therefore  suffice  for  most  engineering  applica¬ 
tions. 


15.  The  recent  acquisition  of  detailed  experimental  data  at 
hypersonic  Mach  numbers  at  both  the  Naval  Ordnance  Laboratory 
and  the  Applied  Physics  Laboratory  (references  c  and  f),both 
with  and  without  steady  state  heat  transfer,  made  it  possible 
to  examine  not  only  the  overall  results  ot  the  theory  but  also 
the  validity  of  the  assumptions  made  and  the  use  of  constants 
drawn  from  incompressible  flow  results. 

16.  The  first  step  is  to  examine  the  conditions  at  the  edge  of 
the  laminar  sublayer.  This  is  necessary  because  the  theory  is 
focused  on  this  point.  For  incompressible  flows  it  has  long 
been  assumed  that  the  value  of  u£  -  y^  at  the  edge  of  the  lam¬ 
inar  sublayer  (in  the  logarithmic  velocity  profile  representation) 
is  roughly  a  constant  that  lien  between  11,0  and  12.0.  The 
present  analysis  is  not  based  on  this  assumption  but  is  so  con¬ 
structed  that  a  computation  and  check  of  the  results  obtained  for 
this  point  may  be  made . 


i 


> 
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17.  The  theoretical  variation  of  uj,  with  Re#  is  compared  with 
experimental  data  (references  i  and  .i)  for  incompressible  flow 
in  Figure  4.  It  is  apparent  that  the  accepted  presumption  that 
u£  -  y£  is  a  constant  is  not  far  from  true;  however,  the  theory 
and  experiment  indicate  that  Reynolds  number  does  influence  the 
value  of  t*L  ~  yt,  slightly. 

18.  Figure  5  shows  a  comparison  between  the  theoretical  and 
experimental  values  (references  c,  f,  and  k)  of  u£  -  y^  for 
compressible  flow.  For  this  comparison  all  values  of  -  yL 

are  based  cn  wall  properties.  The  theoretical  curves  associated 
with  each  set  of  experimental  data  were  calculated  such  that 
they  encompassed  the  experimental  Mach  number  and  Reynolds 
number  range.  Although  the  present  analysis  does  not  accurately 

predict  the  numerical  values  of  -  yi  it  does  predict  the 
proper  trend  of  the  data  for  both  the  heating  and  cooling  cases. 

19.  In  the  formulation  of  this  analysis,  the  incompressible 
skin  friction  law  of  Karman-Schoenherr  has  been  used  as  a  basis 
for  the  zero  heal  U'auuf «r  case.  Unfortunately,  little  exper¬ 
imental  data  are  available  which  describe  the  influence  of  heat 
transfer  on  incompressible  skin  friction  coefficients.  The 
analysis  can,  however,  be  applied  to  this  case  and  comparisons 
made  with  the  fe'w  data  that  are  available.  The  experimental 
results  of  reference  (1),  while  not  reported  in  sufficient  de¬ 
tail  to  make  an  exact  computation  using  the  present  analysis, 
may  be  used  to  show  that  it  does  predict  qualitatively  correct 
results.  Ip  using  these  experimental  data  for  this  comparison 
it  is  assumed  that  the  ratio  Cf/cf<  can  be  used  interchangeably 
with  Cp/C) and  that  little  Reynolds  number  dependence  on  this 
ratio  exists.  Figure  6  shows  both  the  predicted  _:'d  experimental 
variation  of  cf/cy^  with  wall  temperature  ratio  for  a  constant 
value  of  Reg  which  represented  a  mean  for  the  data  of  reference 
(I).  The  results  shown  in  this  figure  demonstrate  that  the  pre¬ 
sent  analysis  describes  coi.vectly  the  variation  of  cf/c#A  with 
increasing  wall  temperature  ratio.  Figure  7  shows  the  influence 
of  heat  transfer  parameter  on  the  values  of  cf/cj^  for  several 
values  of  Mach  number  and  a  single  value  of  Re^.  It  is  ap¬ 
parent  that  cooling  of  the  surface  results  in  an  increase  in 

the  skin  friction  coefficient,  whereas  heating  has  the  opposite 
effect.  That  this  result  is  consistent  with  the  incompressible 
results  shown  in  Figure  6  is  evident  from  the  results  shown  in 
Tigure  8.  This  figure  shows  the  variation  of  cf/c^  with  Mach 
number  for  several  constant  values  of  Tw/T .  The  curves  of 
constant  T^/Tqo  intercept  the  zero  heat  transfer  curvo  at  onlv 
2n<Lpoint‘  JE*0*1  intercept  occurs  at  the  Mach  number  where 
^w/Tcp  "  Te/Too  .  A  curve  of  tho  same  appearance  has  been  de¬ 
duced  by  Seiff  (reference  m)  from  an  empirical  correlation  of 
experimental  data  for  Mach  numbers  up  to  about  5,  A  direct 
comparison  between  the  rosults  presented  in  Figure  8  and  those 
reported  in  reference (m)  is  not  valid  because  the  curves  of 

8  w?r®  computed  for  constant  Rea  and  depend  somewhat  on 
this  Reynold*  number,  whereas  those  ofw  reference (m)  are  as¬ 
sumed  to  l;e  independent  of  Reynolds  number  based  on  distance 
from  the  leading  edge. 
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20.  A  comparison  between  specific  values  of  experimental  skin 
friction  coefficients  and  the  associated  values  predicted  by 
the  present  analysis  is  shown  in  the  following  three  figures. 
Figure  9  shows  the  variation  of  Cf/cf*  with  Mach  number  for 
zero  heat  transfer.  All  of  the  experimental  data  (references 
a,  b,  c,  f ,  k,  n,  o,  and  p)  shown  were  normalized  to  a  constant 
fte$  of  8000  and  are  either  specifically  for  the  case  of  zero 
heat  transfer  or  were  linearly  extrapolated  to  the  zero  heat 
transfer  condition.  The  scheme  used  for  processing  the  heat 
transfer  results  is  outlined  in  reference  (c).  Good  agree¬ 
ment  between  theory  and  experiment  is  found  for  the  entire 
range  of  Mach  numbers  for  which  experimental  data  are  available. 
Figures  10  and  11  show  comparisons  between  theoretical  and  ex¬ 
perimental  values  of  cf/cj ^  plotted  as  a  function  of  heat 
transfer  parameter  for  those  experimental  data  taken  under  con¬ 
ditions  of  steady  state  heat  transfer  (references  c,  f,  and  k). 
For  each  of  the  sets  of  data  shown,  the  variation  of  Reynolds 
number  with  heat  transfer  rate,  if  any,  was  considered  in  the 
theoretical  calculations.  It  is  significant  to  note  that  the 
results  shown  for  the  data  of  reference  (c )  indicate  xxiixe 
variation  of  skin  friction  ratio  with  increasing  heat  transfer. 
From  the  present  analysis  it  appears  that  the  increase  in 
Reynolds  number,  which  accompanied  the  increase  in  heat  trans¬ 
fer  rate,  so  influenced  the  results  as  to  obviate  any  increase 
in  skin  friction  ratio.  In  general,  the  agreement  between 
theory  and  experiment  is  satisfactory  for  each  of  the  sets  of 
data  shown. 


CONCLUDING  REMARKS 

21.  A  theoretical  investigation  of  compressible  turbulent 
boundary  layers  with  heat  transfer  has  been  conducted.  This 
investigation  is  based  on  a  simple  flow  model  which  is  realistic 
for  both  the  zero  heat  transfer  %nd  heat  transfer  conditions. 

The  validity  of  the  flow  model  a>.-.vu*ed  is  demonstrated  by  com¬ 
parisons  between  theoretical  and  experimental  results.  The 
tneory  is  presented  in  such  a  fashion  that  values  of  skin  fric¬ 
tion  may  be  calculated  when  either  the  Reynolds  number  based 
on  boundary  layer  momentum  thickness,  total  boundary  layer 
thickness,  or  distance  from  the  leading  edge  is  given.  Good 
agreement  is  demonstrated  between  theoretical  and  experimental 
values  of  skin  friction  coefficients,  for  both  the  zero  and 
heat  transfer  conditions.  It  is  shown  that  the  predicted 
influence  of  heat  transfer  and  Reynolds  number  on  the  proper¬ 
ties  at  the  edge  of  the  laminar  sublayer  is  consistent  with 
the  available  experimental  data  for  both  incompressible  and 
compressible  flows.  It  is  anticipated  that,  with  the  acquisi¬ 
tion  of  additional  data  covering  a  broader  range  of  conditions, 
improvements  will  be  made  in  both  the  functional  nature  and 
accuracy  of  the  analysis. 
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APPENDIX  A 


General  Relations  for  the  Local 
Skin  Friction  Coefficient 

1.  The  equation  of  motion  to  be  satisfied  by  a  boundary 
layer  flowing  on  a  flat  plate  in  the  absence  of  a  pressure 
gradient  is: 

U  <!h  (A1) 

dx  w  dy  /  dy  J>  dy  dy  '  V  dy/  J 

The  total  shear  stress  at  any  point  in  the  boundary  layer  is 


**!«..  +  *V»b.  -/*£♦/**(&)’ 


Donaldson  (reference  e)  assumed  that  the  ratio  of  the  total 
stress  to  the  laminar  stress  at  the  edge  of  the  laminar  sub¬ 
layer  is  a  constant, 


^Lam.  ^Turb. 
am. 


r  -  constant  at  y  -  S 


Using  this  assumption,  and  evaluating  Eq.  A2  at  y  -  6^  by 
using, 

1  -  n 

du  uod 

tTy  ”  TT7n  (/ 


which  is  obtained  from  the  power  profile  (Eq.  1),  it  can  be 
shown  that  the  thickness  ratio  of  the  laminar  sublayer  is 

n 

2 '  «  *•! 

6l  ^  n(r-l)  SL  (A5) 

6  i2  Re6 

2.  It  may  also  be  logically  assumed  that  at  y  «  Gb,  the 
laminar  shear  stress  is  equal  to  the  turbulent  shear  stress 


ri*«.  ‘  rTurb.  “'-'l 


Evaluating  rTurb<  usiug  Eq.  1  and  TLam  using  Eq.  2  yields  the 
following  expression  for  the  thickness  ratio  of  the  laminar 
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I 

I 

i 

' 

i 

i 

i 

i 


sublayer 


S 


‘n2  5“ 


(A6) 


Equations  (A5)  and  (AG)  necessarily  contain  the  so  far  un¬ 
determined  mixing  length  "j£"  for  which  either  the  Prandtl 
mixing  length  law 


X-  k  y 


or  the  von  Karman  mixing  length  law 


/" 


du 


may  be  used . 


(A7) 


(A8) 


3.  It  is  evident  therefore,  that  the  assumptions  used  to 
obtain  equations  (A5)  and  (A6)  together  with  the  mixing  length 
laws  given  in  equations  (A7)  and  (A8)  will  yield  four  equa¬ 
tions  for  the  thickness  ratio  of  the  laminar  sublayer.  These 
are  tabulated  below; 


Assumption 

_ _ _ _ _ J 

Mixing  length 
law 

- - 

Laminar  i 

sublayer  thick¬ 
ness  ratio  | 

Equation 

(A5) 

Prandtl 

n 

(Eq.  A7 ) 

Sl 

[  n(r-l) 

j 

n  +  1 

6 

L*8  R«*  . 

n 

2n 

von  Karman 
(Eq.  A8) 

SL 

tnri 

J“n(r-l)'j  fl-a 

i  -t-1 

& 

l.*2  «ej 

l  "J 

Equation 

(A6) 

Prandtl 

n 

(Eq.  A7) 

SL 

n  -t-1 

6 

|k*  Re. 

von  Xara&n 
(Kq.  A8) 

n 

2n 

-i  »2  rv* 

j  n  +  1 

.  . — 

6 

|k  2  R*J 

l n 
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The  local  skin  friction  coefficient  for  incompressible  flow 
may  be  calculated  using  the  substitutions 


rw ,  uL 


UL  (&uY 


2  T'w 

cf  -  ± _ „  o 

P<X>  UflD 

together  with  each  of  the  equations  for  the  laminar  sublayer 
thickness  ratio.  These  are  tabulated  below. 


Equation  (A6) 


Assumption 


Equation  (A5)  I  Prandti 
(Eq.  A7 


Local  skin  friction  coefficient 
law 


von  Karmai 
(Eq.  A8 ) 


*25.  JL  2(l-n) 

Cfi  -  2  n<r~*>  ‘  n  JL  ”  +  1fl=B  “ 

L  k2  _Re6.  I  n  (A  10) 


l-o  2 

■iTT-n  r  i  . 


Prandti  f  2  A  i  Tn-t-i 

(Eq.  A7)  Cf ,  -  2  J  1 

1  k7  Ree 


(All) 


*l£  JL  2(l-n) 
von  Karma n  -ilNr  n  1  r  1  .i  n 

(Eq.  A8)  Of.  *»  2  JL  *“n 

[k2  [Re*  pT 


(A  12)| 
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APPENDIX  B 


Derivation  of  Local  Skin  Friction  Coefficient 
Law  lor  Compressible  Turbulent 
Boundary  Layers 

1.  It  has  been  shown  that  equation  (4)  was  the  most  suitable 
for  incompressible  turbulent  boundary  layers.  This  law  is 
Riven  as 

1-n  2 

r  1 1  n  r  in  i-  1 

Cf.  -  9  9. On  I  1  /.o 


To  extend  this  relation  to  compressible  flows  the  property 
values  (density  and  viscosity)  need  be  evaluated  at  the  edge 
of  the  laminar  sublayer.  Equation  (4)  may  be  written  in  the 
form 

1-n  2  _2_ 

„  r  U  +  nr  u  in  i- 1  r . r  i n  >  !|  a  i 


-  2 

20n 

1-n 

I  n 

1 

Ro6 

/go  U0t?  "  0n  [/<*>  Uoo  *j 

using  the  following  substitutions 


A*  ao 

JCI)  Uoo 


in  equation  (Bl),  yields  the  following  relationship: 

1-n  2  i 


-  2 


(20n1 1  f  n  1  n  »  1  p 

L  [Bet>  ; 


n*2U"l 
TU)  nil 

Tl 


2.  To  evaluate  ^oo  it  will  be  assumed  that  the  Crocco 

Tl“ 

quadratic  form  for  the  temperature  distribution  given  as 


T  -  A  |. 


ii 

f-  c 

u 

uoo 

is  valid  By  evaluating  equation  (B3)  with  the  usual  boundary 
conditions,  the  following  equation  results 


I 


NAVORD  Report  3854 


3.  Equation  (B4)  may  also  be  written  as 


(D4) 


with  the  velocity  ratio  at  the  ed^e  of  the  laminar  sublayer 


4.  For  all  calculations  using  equations  (B2),  (B5)  and  (B6) 
or  has  been  taken  as  0.76  and  the  value  of  r.f.  as  0.896. 
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APPENDIX  C 

Derivation  of  Local  and  Average  Skin  Friction 
Coefficient  Laws  Based  on  Re* 

1.  For  incompressible  flow  the  momentum  equation  for  the 
boundary  lay^-r  in  the  absence  of  a  pressure  gradient  is 

d6  „  Pf 

dx  2  (C 


dRee  Pt 

- -sr  (C2  ] 

dRex  2 

The  following  relation  results  when  the  right  hand  side 
of  equation  (4)  is  inserted  into  equation  (C2); 


dReg 


x-n  y,  2 

Hli  n  f  1  n+1  f  n  n  +  1 

- n  +  i  (C3) 

lR  0  l(a  +  l)(n42)J 


and  the  substitution  6 

6 


is  used. 


Equation  (C3)  can  be  integrated  by  assuming  that  n  is 
a  constant  along  a  given  plate  (or  9  is  constant).  Whether 

or  not  this  is  a  valid  assumption  "'ill  be  verified  by  a  com¬ 
parison  between  the  von  Karman  incompressible  mean  skin  fric¬ 
tion  law  and  the  derived  relationship.  Integrating  equation 
(C3 )  yields 

n  f  3  n  -  1  .  2 

n  +  I  ,IL±l[jL]n+1  _ n_ _ T  +  1  Rox 

{  ee;  n  4-1  [20n  (n  ^  i)  (n  +  2) 


a«d  since 


£Fi  _  Re e 
2  Rex 

equation  (C4)  can  be  manipulated  to  yield 
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n  +  1 


A  comparison  between  the  values  of  Cp*  obtained  from 
equation  (C6)  (using  the  n  values  associated  with  the  Ree 

values  of  figure  3  and  equation  (C4))  and  the  values  of  Cpj 
from  the  von  Karman  mean  skin  friction  equation  for  incom¬ 
pressible  flow 

0.242  r_  _  T 

“  LOB10  LCFi  RexJ  CC7) 

indicated  t'^at  over  the  Reynolds  number  (Rex)  range  from 

0.5  x  106  to  60  x  106  the  agreement  between  Cpj  values  is 
within  4  percent.  This  result  means  that  at  least  for 
incompressible  flow  the  assumption  that  n  is  a  constant 
along  a  plate  can  only  lead  to  small  errors  in  the  estima¬ 
tion  of  skin  friction  values. 

Using  this  information,  the  incompressible  law  for 
local  skin  friction  coefficients  on  %  Rex  basis  may  then  be 
derived.  The  momentum  equation  can  be  written  as 


d6  „  6  « 
3S  S3  pt 

Substituting  equation  (4)  yields 
1-n  2 


(C8) 


d«  _  6  L  1  1-n  f  1  In,-1  T  n  ]n+l 

3x  ’  Z5  [2"»J  [(iViHSTzjJ 


(C9) 


which  mav  be  integrated  to  yield  the  foliowing  relation 


°  ±  1  n-1  2 


*n  4  3 

11) 

n  -f  3 

i 

n  f3 

1 

n  -f  3 

(CIO) 

n  -f  1 

■ 

20n 

Rex 

By  inserting  equation  (CIO)  into  equation  (C3)  the  following 
local  skin  friction  law  results 
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•7 

*  .  ■* 


on  Res  or  Rex.  Using  this  assumption,  the  following  relations 
may  be  derived  fro*  an  integration  of  equation  (B2). 


and 


n  i  3 

nTl 


n  + 1 


n  41 


1 

20n 


”  ”  1  fel2  p»r-2-“  r,  f 

1*1  liTJ  N1 


n  +  J 
(C15) 


- 

FTTT 

«4l  n  +  1  n  -  1  n-2 . 52 

m 

Rex 

(»)  („+i)  (20;)  (t^  J 

L  tr 

»  »• 

He»  ■ 


where  Rex  is  given  by 

Rex  -  !  Rex  (maximum) 


n+  3 


(C16) 


To  determine  local  skin  friction  coefficients  with  the 
knowlege  of  only  values  of  Rex  it  is  necessary  to  have  avail¬ 
able  the  variation  of  n  with  Rex  and  Re*  for  use  in  equation 
(C16).  A  calculation  of  the  variation  of  n  with  Rex  and  Re* 
has  been  carried  out  for  the  zero  heat  transfer  case  and  the 
results  are  plotted  in  figures  12  and  13  respective} > .  These 
calculations  vere  made  using  the  results  tabulated  in  Table  I, 
figure  3,  and  equations  (C13)  and  C14).  Curves  similar  to 
those  shown  in  figures  12  and  13  may  be  prepared  for  different 
heat  transfer  rates  by  the  same  method. 


Using  the  data  plotted  in  figures  12  and  13  the  procedure 
for  calculating /Cf  as  a  function  of  Rex  is  as  follows: 

1.  Determine  the  maximum  value  of  Rex  to  be  encounted  for 
a  given  problem.  The  selection  of  this  value  determines  n, 

Re* ,  and  Rex . 

2.  Solve  equation  (C16)  for  T(P  . 

tit 

3.  Insert  the  values  of  u,  and  Re*  in  equation  (B2) 

and  solve  for  Cf.  TL  ' 


By  a  similar  procedure  equation  »C15)  may  be  used  to 
calculate  mean  skin  friction  coefficients. 


As  a  check  on  the  approximate  analysis  described  above, 
the  variation  of  Cf/c^  with  Mach  number  has  been  carried  out 
for  three  values  or  Rex  sufficient  to  encompass  the  Rex  range 
of  the  available  experimental  data.  Figure  1*1  shows  the 


19 


NAYORP  Report  3854 


results  of  these  calculations  compared  to  the  available  ex¬ 
perimental  data.  The  good  agreement  between  theory  and  ex¬ 
periment  tends  to  justify  the  approximations  made  in  the 
preceding  analysis. 

Figure  15  shows  the  variation  of  the  skin  friction  ratio 
with  Mach  number  for  a  number  of  constant  values  of  wall  to 
free  stream  temperature  ratio  and  a  constant  value  of  Rex  of 

107.  The  similarity  between  curves  of  this  figure  and  the 
results  plotted  in  figure  8  is  apparent.  However,  differences 
in  the  curve  shapes  and  corresponding  values  of  the  skin 
friction  ratio  can  be  seen  by  careful  examination. 


20 


V 


NAVORD  REPORT  3854 


TABLE  I 


r~ 

vabiatiom  ur  miHAcnurr  tmcniM  mtio.  4i. 

VARIATION  00  BtRRiACIMWT  TRICtMOO  RATIO, 

T. 

MUUtirttM  nttcRULM  RATIO 

.  1.  ami  MWTOAit  urn 

WWTW  THICK* IRS  RATIO,  |,  A44  WHMUART  UfM  | 

•HAH 

MUNSTMt,  -V"*- 

•  HR  MACN  RlffM*.  *, 

ARB  MAT  HURON* 

•RAM  NJUMtfT 

Ml  91711  RACK  ROMMS,  V, 

A44  4 RAT  HURON*  j 

tOR  A  ViiOCtTV  MoNU  UPOMtRT,  ft, 

,00  1.0 

MRA4074R, 

i,  rot  a  Mucin  monu  unmut,  *,  or  t.o  | 

* 

0 

T.  -T, 

T 

1 

N 

V 

T« 

*.0 

0 

0 

0.10741 

o.uo»i 

1.4000* 

7.0 

• 

0 

0,11*0* 

4.44741 

1.10*74 

t.O 

0.19470 

V. 00140 

3.07313 

1.0 

0.12779 

0.07441 

1. 40741 

4.0 

0. 44011 

0,01*41 

7.09313 

4.0 

0. )••*»• 

0.0*137 

7. 4*1*7 

0.0 

0.40314 

0.03740 

U.  101*1 

0.0 

0. 33*14 

0.03474 

11.40143 

1.0 

0.49441 

0. 01143 

17.34*31 

R.O 

0.03*07 

0. 02414 

*4.44414 

10. 0 

0. 74VOO 

0. 01014 

41.03911 

10.0 

0.71001 

0.017*4 

40.43190 

lt.U 

0.90MB 

0.01314 

*0.40414 

11.0 

0.7*197 

0.01324 

17.34)31 

14.0 

41.43044 

0.01044 

00.300)1 

U.O 

0.400)1 

0.01019 

77.403*0 

10.0 

0.S434T 

0.00470 

104.34744 

U.O 

0,13011 

4.00*91 

101.10040 

u.o 

0.4*173 

0.00671 

131.10413 

U.O 

0.43310 

0.00070 

117. 33930 

10.0 

0.447*7 

0.0011* 

141.74370 

10.0 

0.07130 

0.00140 

IF*  709*3 

1.0 

t.O 

0 

0.11140 

0.00104 

4.04071 

7.0 

1.0 

0 

4. 10010 

0.07114 

3.70470 

t.O 

0.39474 

0.04994 

0.07034 

t.O 

0.33010 

0.00141 

0,34110 

4.0 

0.11041 

0.04937 

10.14041 

4.0 

0.4)370 

0.04*03 

14.01111 

0.0 

0.43941 

0.03370 

14.417*3 

0.0 

0.01*00 

0.0314* 

IT. 9*310 

1.0 

0.70031 

0.03374 

IO.IT744 

*.o 

0.11)44 

0,0114* 

1* .*09 ••  i 

10.0 

0,70737 

0.01731 

44.30710 

u.o 

0.71131 

0.01474 

41.41443 

11.0 

0.00070 

0.01304 

41.01131 

11.0 

0.74711 

0.011*1 

*#, 49190 

14.0 

0.44137 

0.01017 

•1.73313 

10  0 

0.903M 

0.01004 

40.07949 

14.0 

0.04143 

0.00411 

100.70777 

10.0 

o 

0, 00004 

193.47U3 

11.0 

0.40410 

0.00441 

133.7*144 

u.o 

0.0*411 

o.ooo** 

114.4*011 

10.0 

0.40444 

0.00*47 

144.31444 

10.0 

0.1714* 

0.00*40 

lll.OIMt 

«.• 

4.0 

0 

0.41444 

0.M414 

0.00941 

T.O 

4.0 

• 

0.3*411 

0.01701 

0.14091 

t.O 

0.49907 

0.04704 

•.Rill* 

t.O 

0.40004 

9.0*104 

7.74*41 

4.0 

0.44117 

0.04197 

13,00140 

4.0 

0.44414 

0.03471 

11.44111 

4.J 

0.41001 

0.03041 

41.114*1 

4.0 

0.04901 

0.02004 

10.147*1 

0.0 

0.71130 

0.01114 

31.3)000 

0.0 

0.4*417 

0. 01140 

11 .14913 

10.0 

0.77437 

0.014*1 

44.701*4 

to,. 

o.oiiu 

41.1)401  | 

1R.0 

0.41400 

0.01344 

44.40001 

11.0 

0,77100 

0. 01141 

41.11141 

14.0 

0.14417 

0.00991 

0*. 1*404 

14.0 

0.1071* 

0,00070 

41.44179 

1«,0 

0.49731 

0.O0794 

104.13174 

U.O 

0,1340ft 

0.00740 

101.79714 

u.o 

0.01141 

0.00449 

134,4311'* 

u.o 

0.4)414 

0.00040 

111. 14309 

so.o 

0,44070 

0.(0140 

144.439*3 

to.o 

0.07319 

0.00*41 

141.47474 

t.O 

0.0 

0 

0.44447 

0.443*4 

0.11401 

7.0 

0.0 

• 

0.41400 

4.04410 

4.04997 

t.O 

0.11137 

0.044*4 

10.73944 

t.O 

0.41131 

0.04404 

10.13117 

4.0 

0.10341 

0.03114 

3ft. 147*0 

0.0 

0.13*01 

0.03*4* 

14.11009 

0.0 

0.49930 

0.01034 

13.40014 

0.0 

o.ootoi 

0.01404 

31.90039 

1.0 

0.73110 

0.03101 

34.44131 

0.0 

9.47*10 

0. 02  011 

33.10017 

10.0 

3.74047 

0. 01*41 

49.34*14 

10.0 

0.73373 

0.0114* 

47.40107 

11.0 

0.41410 

0.01114 

44,84441 

u.o 

0.77*7* 

0.01103 

04.17100 

14.0 

0. 04140 

0.00947 

IT. *4014 

14.0 

0.R103I 

0.009*7 

14.0*130 

10.0 

o.ooou 

0,00770 

111.14011 

u.o 

0.13479 

0,00770 

107. *3 10* 

u.o 

0.04470 

0.00430 

134.701*0 

1*.4 

0.417*0 

9.000)0 

134.14109 

10.0 

0.00070 

O.OOS33 

1*9. 900*7 

10.0 

0.174*4 

0.90*34 

103,001*9 

1.0 

R.O 

0 

0.03907 

0.04004 

11.03411 

T.O 

1.0 

• 

4.4444# 

0.04*4* 

10.0471* 

t.O 

O.BOIOO 

0.04141 

13.04140 

t.O 

0.413*1 

0.03047 

11.00370 

4,0 

0.43043 

0.03430 

11.03701 

4.0 

0.33413 

0.03230 

17.11*10 

0.0 

0.44141 

0.03430 

10.07*30 

0.0 

0,41471 

9.01*11 

14.01*19 

».o 

0.74100 

0.01940 

37.10014 

R.O 

4.80TS7 

9.019*1 

34.07709 

10. • 

0.7*40* 

0.01*33 

31.70171 

10.0 

0.74000 

t.OUIl 

40.93*49 

le.o 

0.41104 

0.01144 

49.31703 

11.0 

0.71114 

9.01UI 

•O.Otlll 

u.o 

0.04443 

0.00043 

•0.07413 

u.o 

4.41334 

0.00934 

•7.0977* 

u.o 

0.47001 

0.00743 

114.1413* 

u.o 

0.43m 

9.907*1 

110.130*3 

u.o 

O.B47M 

0. 0040* 

141.10037 

u.o 

0.81*14 

0.0001* 

13*. 013*7 

10.0 

o.touo 

0.00)14 

171.409*4 

to.o 

0.47177 

0.00*17 

100,1**90 

1.0 

10.0 

0 

t. O 

4.0 

0.0 

1.0 

10.0 

u. o 

14.0 

U.O 

0,47311 

0,04300 

0,043*1 

0.70001 

0.7100ft 

0.703*0 

0.01174 

0.01100 

0.17141 

0.04044 

0.03774 

0.03111 

0.01411 

0,01444 

0.01443 

0.011*0 

0.00410 

0.0074# 

14.11410 

13.73003 

>0.03041 

10.1*077 

30.7711* 

14.1703* 

71.00000 

oi.ioru 
no. ootoo 

141. 0*109 

T.t 

10.0 

• 

t.O 

4.0 

4.0 

1.0 

10.0 

u.o 

14.0 

U.O 

u.o 

•  .Mill 
4.4171* 
4,440(4 
4.44143 
4.44931 
0.74704 
0.711)1 
4,11411 
4.44000 
0.14073 

4. 4370* 
4.43*44 
4.00471 

0,413*3 

0.01*31 

0.0U3I 

0.0113* 

0.0001) 

0,00747 

0,00410 

19.11419 

14. *7009 

10. 031)9 
17.30131 
30.174*9 
*t.U4*0 
99.10700 
00.40**4 

111. *702* 

13*. *97*0 

to.o 

0.00141 

0.00*19 

173.09*10 

to.o 

0.8741) 

0.00:1 

i 


NAVORD  REPORT  3854 


TABLE  I  (CONTINUED) 


VARIATIOR  or  OIIPUCmiTT  TNICKMM  RATIO, 

•OAWTtW  TIUCriMUt  MTIO,  |,  ADO  *099  MAT  UTH 

iwn  nuunrra  wit*  iuoi  wumica,  m,  add  niat  TMiim 

MMttkTM,  T«t^T»,  ro«  A  VfVOCtTT  POOTIL*  *1*09*147,  «,  t)»  ».• 

tui, tiw  or  ouruciauT  meant  uno,  41, 

■Do»rw  mean,  uno,  |,  w  Kow.t  uno 
urn  Nuna  tin  out  ttoon,  wo  out  nunm 

MUOtTU,  TOO  A  T.L0C1TT  fOOTIU  UfoWT,  a.  Of  11.. 

9 

* 

*• 

*m 

■£ 

• 

1 

n 

0 

T.  -  T. 

* 

1 

N 

"T.  ' 

».o 

0 

0 

0.10001 

0.00  loo 

1.1*071 

11,9 

9 

9 

0.99171 

9.9704* 

1.1IM* 

M 

0.10041 

3.0*003 

*.7*077 

t.9 

9.103*9 

t. 73373 

4.0 

0. 31100 

0.04714 

7.4*00* 

4.0 

9.3146* 

7.93333 

0.0 

0,41001 

0.03331 

10.1**40 

0.9 

9.41030 

0.03011 

14. *0311 

B.O 

0.10314 

V.OIIOO 

a*. *io*i 

9.9 

O.IIUf 

10. « 

0.00011 

0.01*00 

30.0*000 

U.O 

0.03437 

11.0 

0.71164 

0.01*01 

i*.4f *14 

ia.o 

0.0*41* 

0.01*44 

14.0 

0.70030 

0.01007 

70.100*3 

14.0 

0,74000 

1«.0 

0.00131 

0.00000 

•*.1707* 

19.0 

0.77104 

10.0 

0.0III4 

0. 00001 

114, *14*0 

19.9 

0,00410 

0.00000 

III. 7130* 

10.0 

0.14700 

0.00111 

1*3. 1*3*1 

*9.9 

0.410*0 

0.00141 

161,904*9 

1.0 

1.0 

0 

0.10170 

0.00103 

s.tout 

11.9 

t.9 

9 

O.lOftll 

0.03*77 

3.90199 

S.O 

o.moo 

0.0*107 

1.1311* 

1.0 

9.94*31 

0.04*01 

9.0*013 

4.0 

0,40177 

*.04117 

• .7*0*1 

4.0 

0.303*0 

0.0371* 

«.o 

0.11710 

0. 0*1*1 

IT. 4, TOT 

9.0 

0.4773* 

0.03777 

17. 1**4* 

t.o 

0.00011 

o.oauo 

1*. 1*741 

9.0 

0.37105 

10.0 

0.07IS1 

0.01*17 

41 ,*41*7 

10  0 

0.043*1 

0.01314 

11.0 

0,73117 

0.01140 

U.73BM 

11.0 

0,70031 

0.01107 

90.01*1* 

14.0 

0.7TI07 

0.00*01 

31.14117 

14.0 

0.74411 

0,00*3* 

10.0 

i.iun 

0.007*1 

101.10**3 

19.0 

0,77(04 

11.0 

9.11909 

o. 00*11 

117 ,1*10* 

11.0 

0.00014 

so.o 

0.04*11 

0.00144 

110.10**4 

*0.0 

0.03630 

0.00340 

1)3.3111* 

1.0 

4.0 

• 

O.lttl" 

0.011*4 

l.***4* 

11.9 

4.0 

0 

0.17171 

0.04700 

». 7*0*9  S 

1.0 

0.11034 

*. 04*01 

.7,47701 

9.9 

0.31147 

0.04*7* 

7.30311 

4.0 

0.44177 

0.03171 

It.  com 

4.9 

0. 40301 

0.0 

0.04077 

0.01730 

11.7*000 

9.9 

0.9013* 

19.43411 

0.0 

0.01110 

1.0SO4* 

30.47*** 

9.9 

0,0*330 

0.01*4* 

30.036*1 

10.0 

0.00070 

0.01*03 

44.117*4 

10.0 

0.06313 

**  **  1 4  *0 

ia.o 

0,73070 

0.01*07 

*1.0310* 

lt.0 

0.70004 

<0. 93**1 

14.0 

0.773*0 

0.00000 

lo.aotoo 

14.0 

0. 74101 

0.00*3* 

7*. *1433 

10.0 

0,00010 

0.00770 

103.0003* 

19.0 

0.71134 

11.0 

o.oioot 

0.00041 

i»*.  ooom 

J«.0 

0.00*00 

SO.O 

0,0104* 

0. 00037 

11*. 37*1* 

Zv.O 

O.M3070 

0,00)31 

190.14310 

M 

l.l 

0 

O.SOOIf 

0.044*1 

I.IMlft 

11.0 

9.9 

9 

9.31*43 

0,040*1 

S.043H 

1.0 

0,40117 

0.0410* 

*.7*147 

t.9 

0.3*1*0 

0.037*4 

4.0 

0.47771 

0.03311 

i, .si.it 

4.0 

0.43736 

0,03101 

0.0 

0.10141 

0.0*147 

*1.04114 

9.0 

0,13341 

0.0*410 

fl. 0/039 

f.O 

0.03100 

0.0111* 

31.74**4 

9.9 

0, 6**40 

0.01*34 

33.9*159 

0.134*3 

0.014*4 

4*. 40**4 

10.0 

0.0*11* 

0.01441 

46.0*4*0 

11.0 

0.741*1 

0.01173 

0),*4S0t 

lt.0 

0,71166 

0.0113* 

0.77*17 

0.00037 

*3.16*0* 

14.0 

0,70171 

0.00*1* 

9.9M79 

0.007*3 

101.7*14* 

19.0 

0.703*0 

10,0 

0.03314 

0.00031 

131 .*3344 

14.0 

0.100*4 

0.00023 

0.1*17* 

0,0011* 

101.0*07* 

10.0 

0,13107 

0.00533 

19* .01013 

t.o 

1.0 

0 

0.4104* 

0. 03*37 

19.99971 

11.0 

1.0 

0 

0.37340 

t.o 

0.44*1* 

0.03*00 

11.0*74* 

t.o 

0.40306 

O.KOOIO 

o.towr 

10. *77*4 

4.0 

0.44*33 

0.37*** 

0.03314 

*4.31404 

9,0 

0.64133 

0.040*0 

0.01044 

30.03700 

9.0 

0.41013 

0.01770 

0,70*07 

0. 01440 

40.710*4 

lo.o 

0  44004 

11. 0 

0.74*70 

0.01130 

09.IM2* 

13  U 

0.71040 

14,0 

11,70311 

0.00017 

00.331*1 

14.0 

0.7153* 

10.0 

0.01113 

0,00710 

104.14**7 

1B.0 

0, 2**10 

11.0 

0  MJ4I1 

0.00433 

134.11070 

14.0 

0  MUM 

0.1131* 

0.00333 

113. 400*1 

10.0 

O.M334I 

0. 00117 

tfl ,005*0 

*.o 

10.0 

0 

0. 4*301 

0.03134 

It. *43*0 

11.0 

10.0 

0 

0,41330 

0,033*7 

0.47031 

0.033*3 

14.37104 

3.0 

0,431*1 

4,0 

0.31111 

0,07*03 

10.93444 

4.0 

0.4*137 

0,024.43 

4.0 

0,39031 

0.02243 

IB. 40113 

C.O 

0.51HIA 

«  o 

V.tlliM 

O.UI/feJ 

37.34304 

4.0 

u.07ir.n 

0,'ilMl 

10,0 

0.70*07 

0.013*0 

01. 01434 

10.0 

I)  A7MA 

12  0 

0.7014* 

0  niiuH 

*7.13491 

13,0 

0.711*3 

0  01 010 

It  0 

0.71373 

0  0*J«»7 

*7.391*2 

14.0 

0.71*94 

U  o 

0.01340 

0  00734 

110.31*30 

1C.0 

<1,74*1*1} 

0  OU774 

IB  0 

0.BJ5H 

0,(«»tl  1 

13A  34313 

14.0 

0  41310 

0  00104 

0.01417 

O.W317 

141,37344 

30  n 

0  1.3  17  J 

V  4m»>12 

i 

s 

I 


i 

i 

i 

1 


1 

4 

* 


t 

i 


t 

« 

<4 

•4 

•i 

t 

■4 


1 


22 


< 


NAVORD  REPORT  3854 


TABLE  I  (CONTINUED) 


yam  at  io*  or  uu  flat  hunt  tmcikma  iutio,  -J-, 

IMIANTUM  IlllCIOkJI  AATIU,  {,  A6U  DOW  DA  If  UTSJt 

IKAM  rAMIU7U  WITH  HACK  9WUKR.  M,  A6U  NUT  TOOIIU 

T-  -  T» 

NAtANCTU,  «-*,  ft*  A  TilAJClTY  PDOTIU  9X90*1077,  A,  9.0 

VA9U1I0N  Of  DtirtACMm  tnicmu*  iaiiu,  -f, 

■MUTIN  TMCVKUA  RATIO,  J,  A«0  MKTUAAT  LAYER 

IIUK  MUItTO  »1TN  MAfM  KUNti**,  M,  AND  NUT  TRARIfftN 

Nuinn,  T%f'  Tt,  nx  a  yuajcitt  mom*  kyrmutt,  •«  or  7,0 

A 

% 

~T£ 

4-*. 

» ' 

• 

i 

M 

A 

- 

* 

6 

J 

a 

9.0 

-2.0 

4.0 

0.3*103 

0.07799 

5.3*766 

7.0 

-2.0  <4.0 

0.31171 

0.0004 

9.2141ft 

6.0 

0.572U 

0.04219 

13.57533 

6.0 

0.50724 

0.03166 

13.11272 

6  0 

o  r*>ow# 

0.02735 

24.6*5x0 

6.0 

0.4225* 

0.025*6 

*4.0571* 

10,0 

0.71152 

0.0100V 

36.37140 

10.0 

0.700*7 

0.01141 

36.07170 

12.0 

0.60032 

0.01404 

37.00142 

12.0 

0.75632 

0.01372 

55,133*1 

14.0 

0  6V.37 

0.01074 

77.763** 

14.0 

0,70702 

0.01057 

75.4020* 

u.o 

0.6*144 

<M*J*47 

lot. 70017 

16.0 

0.62773 

0.00639 

*6.65316 

11  0 

0.6N13S 

0.00664 

126.63*69 

U.O 

0.65145 

0.0CG61 

135,03671 

20.0 

(1.60666 

(i. DU  503 

159.30726 

30.0 

0,97011 

0.00566 

1*3.73762 

9.0 

•4  It 

i.U 

U.633JH 

0.04643 

11.01363 

7.0 

•4.0  6.0 

0.46*27 

ft  ft4176 

10,70400 

6.0 

0 . 66316 

II.  02*66 

21  376*4 

6.0 

0.60471 

ft.li2T64 

71  . 111  Ml 

lll.U 

U. 7 4263 

0.02013 

36.6*021 

lo.o 

0.49094 

0.01*3* 

39.76097 

lY  U 

11.70919 

0,01457 

54.97790 

12.0 

0.75044 

0.01410 

53.46231 

14.0 

0.63210 

0.01105 

73.31222 

14.0 

0.7*330 

0.01065 

73.11521 

16.0 

0.65034 

0.00665 

»V. 34104 

U.O 

0.92528 

0.00659 

*6.07663 

U.O 

0. 676*2 

0.000*5 

126. 4.0432 

16.0 

0.64974 

0.000*3 

17i.6|U3 

20. U 

0.6*567 

0.00571 

15*. *0016 

20.0 

O.66HV0 

0,00573 

151.90599 

9  0 

n.o 

0.46265 

0.01736 

6.41234 

7.0 

•6.0  6.0 

0.416*2 

0 .05000 

9.26663 

h.O 

O.M36U 

U. 03242 

1*. 65611 

6.0 

»  56376 

(1.0301* 

1*. 370*7 

lu.ft 

0.73260 

0.02131 

34.36701 

10.0 

•».  «•*•»* 

'.02036 

33.40373 

12.0 

0.76*63 

0.01510 

52.06443 

12.0 

0.744U 

0.01473 

50.92274 

14.  * 

0.62674 

0.01137 

72.6*356 

14.0 

0.76*40 

0.01118 

70.73477 

10. .1 

0. 65716 

U.  (>066  4 

*6.96433 

1C.0 

0.92270 

0  00675 

*4.022*6 

U.ft 

(l  67644 

0,00706 

124.06760 

1N.0 

0.6479* 

0.00705 

12(1.2636* 

20.0 

0.64463 

0.0(l56« 

154.37360 

20.0 

0.9C706 

0. 00170 

14*. 9(1*3 

1.0 

•  *.U 

h.O 

U. 62047 

0.01563 

17.3176* 

7.0 

*6  0  6. ft 

0.55*00 

0.03303 

U, *2401 

10.0 

U. 72106 

0.02263 

31. *0433 

10.0 

ft  6C624 

0.02153 

31.0357* 

12.0 

0.76372 

0.01560 

49.60127 

12.0 

V.i.lJ 

0.01531 

46.17113 

U.O 

0.62525 

0.01172 

70.4160V 

14.0 

0.71532 

U. 01146 

69.405V? 

ir,.o 

0.X54KD 

0.00603 

04.46409 

U.O 

0.62005 

0.006*4 

•it  .7337* 

U.O 

0.670*1 

0.00771 

121. *2275 

16.0 

0,6441* 

0.00717 

lU.nU53 

20. U 

O.NV377 

0.00566 

152.00060 

20.0 

0.6*61* 

0.00567 

147.5*7*6 

1.0 

-10.  u 

6.0 

0.50244 

0.0401* 

14,739*9 

7.0 

-10.0  6.0 

0.53041 

0.03656 

14.5076C 

10.0 

0.700*1 

0.024U 

*3.36326 

10,0 

0.63516 

0.02263 

2* .*73*6 

12,0 

0.77736 

0.01650 

*  -.1515 

12,0 

0.73040 

0.015*3 

4ft. *5060 

U.O 

0.62192 

0.0120* 

6,  .  S*72 

14.0 

0,76104 

0.01  IMS 

66.0744ft 

U.O 

0,65253 

0.00*27 

Vl.vYJ32 

u.o 

0.61779 

0.00*15 

69.37240 

U.O 

0.67534 

0.00733 

110.06644 

U.O 

0.64433 

0.00T2B 

115. 6161* 

20.0 

0.6*206 

0,00507 

14V. 5309* 

20.0 

0. 6(5509 

0.005*6 

145. 1  ft  1 01 

23 


\ 


i 


NAVORD  REPORT  3854 


TABLE  I  (CONCLUDED) 


VARIATION  OF  DtSKACm«T  THtCIOtm  RATIO,  -J-, 
wmrnrn  THicukii  ratio,  |,  and  huvroarv  utar 

RNA n  FARAMPTFR  «ITH  MACH  HUHhCR,  N,  ARM  HUT  TRARIFR 
IAMKTRR,  TmJLJx,  ro«  A  VELOCITY  7ROF1U  At  FOR  CRT,  n,  OJF  1 


VARIATION  Of  BISPUCItniT  TMICVNRRS  RATIO,  4*. 

■ORUrrW  THICRHWS  RATIO,  J,  ARU  KWOARY  UTIR 
•MAH  FMAIRTUt  VITM  UC1I  HIAIPRR.  M,  ARti  DRAT  T4AMFA4 
FARAMiTUt,  Ft*  A  t  LUC  ITT  MOFIU  IlHHUrT,  »,  OF  11. Q 


ll. 0 

•1,0 

4.0 

0.2SIS0 

0.04460 

3.04412 

C..0 

0.4U7U 

0.03294 

12.49172 

R.O 

0.33MO 

0.07300 

23.33044 

10.0 

0.42414 

0. 01449 

36.93046 

12.0 

0.66774 

0.01244 

33.33419 

14.0 

0.73543 

0.01000 

73.14113 

14.0 

0.7)240 

G.OORbC 

93.90371 

U.o 

0.H0217 

0.00060 

121.33030 

20.0 

O.R233I 

0.0*630 

*130.07273 

tl.O 

-4.0 

4.0 

0.3KIU 

0.03634 

10.42100 

R.O 

O.SU74 

0.02439 

21.09344 

10.0 

0.41313 

0.01763 

34.73634 

12.0 

0.69U97 

0.01323 

31.39623 

14.0 

0.73143 

O.OlOJl 

70.94041 

14.0 

0.70994 

0.00423 

91.64013 

U.O 

D.R0003 

o.oor.n 

119.22304 

30.0 

0.42397 

0.00337 

147.91741 

11.0 

-4.0 

0.0 

0.33S4K 

0.04123 

6.13724 

n  n 

0.49707 

0.02639 

U.I364I 

10.0 

0.40131 

0. OU30 

32.49730 

12.0 

0.1*7391 

0.0137U 

40.1424* 

14.0 

0.72443 

0.01037 

41.77010 

U.O 

0.7644? 

n  DOfci- 

91.31351 

u.o 

0.707KV 

0  '<)>••» 

■•6.99413 

20.0 

O.H2233 

0.00364 

I45.797H7 

11.0 

-R.O 

R.O 

0,47244 

0.02130 

14.37193 

lu.o 

0.39921 

0.01944 

30.21133 

12.0 

0.04470 

0.01421 

46.RI243 

14. n 

0.73205 

0.01043 

46.63300 

U.O 

0.74371 

0. OCA 34 

44.21719 

u.o 

0.7934* 

0.00493 

114.41962 

20.0 

O.R2U7S 

w.wv372 

143.49330 

ll.O 

-10.0 

R.O 

0.44414 

0.03104 

14.30733 

lo.o 

0.&741T 

0.02049 

24.02343 

12.0 

0.63411 

0.C1473 

44.67733 

14.0 

0.71702 

0. 01114 

64.30493 

11.0 

0.76040 

0.00474 

R1.0U139 

U.O 

0.79J41 

0,00704 

112.69944 

20.0 

0.41914 

0.003R1 

140.94107 

o?K 


NAVORD  REPORT  3654 


0  20  40  60  60  100 

PERCENTAGE  JF  TOTAL  SHEAR  STRESS 


FIG.  2  SCHEMATIC  VIEW  OF  VISCOUS  AND  TURBULENT 
SHEAR  STRESS  ACROSS  A  TYPICAL 
TURBULENT  BOUNDARY  LAYER 


VARIATION  OF  TURBULENT  BOUNOARY  LAYER  VELOCITY  PROFILE  EXPONENT  (n) 

WITH  REYNOLDS  NUMBER  (Ree) 
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FIG. 5  COMPARISON  BETWEEN  THEORETICAL  AND 
EXPERIMENTAL  VALUES  OF  u+  =  yL+ 

FOR  COMPRESSIBLE  FLOW 
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FIG.  6  INFLUENCE  OF  HEAT  TRANSFER  ON  SKIN 
FRICTION  RATIO  FOR  INCOMPRESSIBLE  FLOW 


FIG. 7  INFLUENCE  OF  HEAT  TRANSFER  ON  SKIN 
FRICTION  RATIO  FOR  THREE  VALUES 
OF  MACH  NUMBER  AND  Ree  *13,500 


FIG.  8  VARIATION  OF  SKIN  FRICTION  fiATIO  WITH  MACH  NUMBER  FOR  SEVERAL 
CONSTANT  VALUES  OF  WALL  TO  FREE  STREAM  TEMPERATURE  RATIO  ANO  Re»l3,500 
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FIG.  10  COMPARISON  BETWEEN  THEORETICAL  AND 
EXPERIMENTAL  VALUES  OF  SKIN  FRICTION  RATIO  FOR 
M«2.43,5,0  AND  6,8 
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FIG.  II  COMPARISON  BETWEEN  THEORETICAL  AND 
EXPERIMENTAL  VALUES  OF  SKIN  FRICTION  RA'HO 
FOR  M»5.75,8.25,  AND  9.0 
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Fi  .13  VARIATIC..  n  WITH  Re*  FOR  SEVERAL  CONSTANT  VALUES  OF  MACH  NUMBER 


REF. 


FOR  CONSTANT  VALUES  OF  Rex  OF  I06, 107,  AND  I08 
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